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Abstract. 



We have reconstructed the three-dimensional density fluctuation maps to z ~ 1.5 using the distribution of galaxies 
observed in the VVDS-Deep survey. We use this overdensity field to measure the evolution of the probability 
distribution function and its lower-order moments over the redshift interval 0.7< z<1.5. We apply a self-consistent 
reconstruction scheme which includes a complete non-linear description of galaxy biasing and which has been 
throughly tested on realistic mock samples. We find that the variance and skewness of the galaxy distribution 
evolve over this redshift interval in a way that is remarkably consistent with predictions of first- and second-order 
perturbation theory. This finding confirms the standard gravitational instability paradigm over nearly 9 Gyrs of 
cosmic time and demonstrates the importance of accounting for the non-linear component of galaxy biasing to 
consistently reproduce the higher-order moments of the galaxy distribution and their evolution. 

Key words, cosmology: large scale structure of the Universe — cosmology: theory — galaxies:statistics — galaxies :high- 
redshift — galaxies: evolution 



1. Introduction 



According to Thomas Wright mapping the Cosmos on 
the very largest scales is about gaining "a partial View 
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of Immensity, or without much Impropriety perhaps, a fi- 
nite View of Infinity '0. Unfortunately, charting the cos- 
mic territory beyond our local volume into the distant 
Universe is observationally challenging. Until recently, our 
understanding of the large-scale organisation of galaxies at 
z > 0.2 had to rely on the predictions of numerical simu- 
lations in the frame work of the rather suc cessful cold dark 
matter model (e.g., ISpringel et al. . 2005 ). 

Within this scenario, which has now developed into the 
leading theoretical paradigm for the formation of struc- 
tures in the Universe, structures grow from weak, dark- 
matter density fluctuations present in the otherwise homo- 
geneous and rapidly expanding early universe. The stan- 
dard version of the model incorporates the assumption 
that this primordial, Gaussian-distributed fluctuations are 
amplified by gravity, eventually turning into the rich struc- 
tures we see today. 

This picture in which gravity, as described by general 
relativity, is the engine driving cosmic growth is gener- 
ally referred to as the gravitational instability paradigm 
(GIP). However plausible it may seem, it is important to 
test its validity. In the local universe the GIP paradigm 
has been shown to make sense of a vast amount of indepen- 
dent observations on different spati al scales from galaxies 



to superclusters of g alaxies (e.g., iPeacock et all 12001 



Tegmark et all l2006f) . Deep redshift surveys now allow us 



to test whether the predictions of this assumption are also 
valid at earlier epochs. 

In this paper we test the role of gravity in shap- 
ing density inhomogeneities by using three-dimensional 
maps of the distribution of visible matter revealed by the 
VIMOS-VLT Deep Survey over the large redshift baseline 
< z < 1.5 (see Massey et al. 2007 for three dimen- 
sional cartography of mass overdensities in the COSMOS 
field). We present first a qualitative picture of the large- 
scale organization of remote cosmic structures, and then 
quantify the observed clustering by computing the proba- 
bility distribution functions (PDF) of galaxy overdensities 
5 g . In this way, we trace how the amplitude and spatial 
arrangement of galaxy fluctuations changes with cosmic 
time. We explore the mechanisms governing this growth 
by comparing the time evolution of the low-order moments 
of the galaxy PDF, (i.e. the variance amplitude < Sg > 
and the normalised skewness S3 —< Sg > c / < Sg > 2 ) 
with the corresponding quantity theoretically predicted 
for matter fluctuations in the linear and semi-linear per- 
turbative regime. (Note that in the following we shall of- 
ten speak equivalently of the variance or of its square root, 
i.e. the root mean square amplitude < Sg >V 2 when refer- 
ring to the second-order moment). This provides a test of 
GIP-specific predictions at as-yet unexplored epochs that 
are intermediate between the present era and the time of 
decoupling. Knowledge of the precise growth history of 
density inhomogeneities pr ovides also a way to test the 
theory of gravitation (e.g., Linder . 20051 ). 



In addition to the statistical approach presented in 
this paper, we have recently addressed this same issue 
also from a dynamical point of view. We have used lin- 
ear redshift-space distortions in the WDS- Wide data 
to meas ure the growth ra te of matter fluctuations at 



0.8 ( Guzzo et all 120081 ) . This approach offers promis- 



ing prospects for deter mining the cau se of cosmic accelera- 
tion in the near future (lLinderl . [2007h . The work presented 
here is also complemented by a parallel paper (Cappi et 
al. 2008) in which we study the behavior of the N-point 
correlation functions for this same sample. Higher-order 
galaxy correlation functions are known to display a hier- 
archical scaling as a function of the variance of the count 
distribution (e.g., Peebles 1980). In the same spirit, we 
use this scaling to test the standard assumption of evolu- 
tion under gravitational instability of an initially Gaussian 
distribution of density fluctuations. 

The paper is organised as follows: in §2 we briefly de- 
scribe the first-epoch VVDS data sample. In §3 we present 
3D overdensity maps from the galaxy distribution in the 
VVDS to z ~ 1.5; we then characterise the evolution 
of galaxy fluctuations with cosmic epoch by computing 
their PDF in two redshift slices. In §4 we compare the 
observed redshift evolution of the low-order moments (i.e. 
variance and skewness) of the PDF of the galaxy fluctu- 
ations with linear and semi-linear theoretical predictions 
of the Gravitational Instability Paradigm. Conclusions are 
presented in §5. 

The coherent cosmological picture emerging from inde- 
pendent observations and analysis motivates us to present 
our results in the context of a ACDM cosmological model 
with Q m = 0.3 and J7a = 0.7. Throughout, the Hubble 
constant is parameterised via h = i?o/100. All magni- 
tudes in this paper are in the AB system (Oke & Gunn 
1983), and from now on we will drop the suffix AB. 



2. The First-Epoch WDS-Deep Redshift Sample 

The primary observational goal of the VIMOS-VLT 
Redshift Survey as well as the survey strategy and first- 
epoch observations in the VVDS-0226-04 field (from now 
on sim ply VVDS-02h) are presented by iLe Fevre et al 
(2005). Here it is enough to stress that, in order to min- 



1 An Original Theory of the Universe (1750, 9th letter, Plate 
XXXI). 



imise selection biases, the VVDS-Deep survey has been 
conceived as a purely flux- limited (17.5 < / < 24) survey, 
i.e., no target pre-selection according to colors or com- 
pactness is used. Stars and QSOs have been a-posteriori 
removed from the final redshift sample. Photometric data 
in this field are complete and free from surface bright- 
ness sele ction effects, down to the limiting magnitude 
Jab =24 (|Mc Cracken et alll2003[ ). Spectroscopic observa- 
tions were carried out using the VIMOS multi-object spec- 
trograph using one arcsecond wide slits and the LRRed 
grism which covers the spectral range 5500 < A (A) < 
9400 with an effective spectral resolution R ~ 227 at 
A = 7500A. The rms accuracy in the redshift measure- 
ments is ^275 km/s. Details on the observations and 
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data reduction are gi ven in iLe Fevre et al.l (|2004f ) and in 
Le Fevre et ail (|2005l) . 

The VVDS-02h data sample extends over an area of 
0.7x0.7 sq.deg (which was targeted according to a one, 
two or four passes strategy, i.e., giving to any single galaxy 
in the field one, two or four chances to be t arget ed by 
VIMOS masks (see fig. 12 in ILe Fevre et al.l (|2005f) l has 
a median depth of about z~0.76. It contains 6582 galax- 
ies with secure red shifts (i.e., reds h ift de termined with a 
quality flag>2 (see ILe Fevre et all (|2005l ))) and probes a 
comoving volume (up to z = 1.5) of nearly 1.5 • 10 6 ft.~ 3 
Mpc 3 . This volume has transverse dimensions ~ 37 x 37 
/i _1 Mpc at z = 1.5 and extends over a comoving length 
of 3060 h~ 1 Mpc in the radial direction. 

For the statistical analysis presented in this paper, we 
first define a sub-sample (VVDS-02h-4) including galax- 
ies with redshift z<1.5 and over the sky region (0.4x0.4 
deg 2 ) that was repeatedly covered by four independent 
VIMOS observations in each point. Even if measured red- 
shifts in the VVDS reach up to z~5 and cover a wider 
area, these conservative limits bracket the range where we 
can sample in a denser way the underlying galaxy distri- 
bution and, thus, minimise biases in the reconstruction of 
the density field (see the analysis in §4.1). The VVDS- 
02h-4 subsample contains 3448 galaxies with secure red- 
shift (3001 with 0.4 < z < 1.5), probes one-third of the 
total VVDS-02h volume and it is characterised by a red- 
shift sampling rate of ~ 30% (i.e. on average about one 
over three galaxies with magnitude Tab <24 has a mea- 
sured redshift). This high spatial sampling rate is a crit- 
ical factor for minimising biases in the reconstruction of 
the 3D density field of galaxies. To optimise the analy- 
sis of the associated probability density function, we fur- 
ther select only galaxies with absolute blue magnitude 
Mb < — 20 + \ogh. With this selection, we define two 
nearly volume-limited sub-samples in the redshift ranges 
0.7 < z < 1.1 and 1.1 < z < 1.5 respectively. A discus- 
sion of possible effects of galaxy evolution on our results 
is presented in § 4.3. 

3. The galaxy density field at high redshift 

The first large redshi ft surv e ys of the local U n iverse 



Davis fc Huchra. 1981 



Giovanelli k. Havnesl 1991 



Geller & Huchra 
Strauss et al.l 



1991 



1992a: 



da Costa et all 1994 ) showed that galaxies have a 



highly non-random spatial distribution and cluster in a 
hierarchical fashion. The corresponding three-dimensional 
maps reveal a complex web-like network of thin, filamen- 
tary structures connecting centrally condensed clusters of 
galaxies, punctuated by large, quasi-spherical, low-density 
voids. These structures are the outcome of more than 13 
billion years of evolution of small-amplitude fluctuations 
that we see reflected in the temperature anisotropy of 
the Cosmic M i crowa ve Background (CMB) at z ~ 1100 



( Spergel et all 120071 ). Recent analyses (e.g. Tegmark et 



al. 2006) have shown the remarkable consistency between 
two-point statistics of the galaxy distribution at z ~ 



and the CMB power spectrum which probes matter 
clustering at the recombination. Mapping the large-scale 
structure at z ~ 1 is thus crucial to further test the 
coherency of the gravitational instability picture at an 
intermediate time between the epoch of last scattering 
and today. 

In this section we present a reconstruction of the 3D 
galaxy density field, discussing first the methodology and 
summarising the techniques adopted to correct for ob- 
servational selection effects. These are fully presented in 
Marinoni et al. (2005, hereafter Paper I) and Cucciati et 
al. (2006), to which the reader is referred for more details. 

3.1. Density reconstruction method 

The continuous galaxy density fluctuation field 



<S fl (r,R) 



p(r,R)-p 



(1) 



represents the adimensional excess/deficit of galaxies on a 
scale R, at any given comoving position r with respect to 
the mean density p. As suggested by Strauss and Willick 
(1995) we estimate the smoothed number density of galax- 
ies brighter than M c on a scale R, p(r,R,< M c ), by 
summing over an appropriately weighted convolution of 
Dirac-delta functions with a normalised Gaussian filter F 



p(r, R, < M c ) 



E 



F(u) = (2irR 2 ) 3/2 exp 



f™6 D (u-\An\/R)F(u)du 



S(ri,M c )®(m)t(ri,m)V(a, 5) 



(2) 



(3) 



Here Ar = (r\ — r) is the separation between galaxy po- 
sitions and the location r where the density field is evalu- 
ated. We compute the characteristic mean density at po- 
sition r using equation (|2|) by simply averaging the galaxy 
distribution in survey slices r±i? s , with R s = 400/i~ 1 Mpc. 
The four functions in the denominator of equation [2] cor- 
rect for various observational characteristics: 

- S(ri,M. c ) is the distance-dependent selection func- 
tion of the sample. This function is identically one when a 
volume-limited sample is used. When the full magnitude- 
limited survey (17.5 < I < 24 in our case) is used, how- 
ever, this function corrects for the progressive radial in- 
completeness due to the fact at any given redshift we 
can only observe galaxies in a varying absolute magni- 
tude range. While the PDF of galaxy fluctuations will be 
derived from volume-limited samples, in the next section 
we shall make use of this function when reconstructing a 
minimum-variance 3D density map from the full VVDS 
survey. 

The actual values of S(r, Ai c ) are derived using the 
VVDS galaxy luminosity function (Ilbert et al. 2005), as- 
suming a minimum absolute magnitude A4 C = —15 + 
5 log h and accounting for its evolution as measured from 
the VVDS itself. A more detailed discussion of the deriva- 
tion of the selection function can be found in Paper I 
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- $(m) corrects for the slight bias against bright 
objects introduced by the slit positioning tool 
VMMPS/SPOC (Bottini et al. 2005). 

- C( r z; m ) is the correction for the varying spectroscopic 
success rate as a function of the apparent Iab magnitude 
and of the distance of the object (see Ilbert et al. 2005). 

- \I/(a,(5) is the angular selection function correcting 
for the uneven spectroscopic sampling of the VVDS on 
the sky (see Fig 1. of Cucciati et al. 2006). Its purpose is 
to make allowance for the different number of passes done 
by the VIMOS spectrograph in different sky regions (a 
factor which is anyway maximised in the 4-pass sub-area 
of the sample) . 

The analytical form of these selection functions is dis- 
cussed in Cucciati et al. (2006). The underlying assump- 
tion in this reconstruction scheme is that the subset of 
observed galaxies (e.g. in the case of a flux-limited sam- 
ple, those luminous enough to enter the sample at a given 
redshift) is representative of the full population. This as- 
sumption clearly neglects any dependence of clustering on 
luminosity and could bias the density field reconstructed 
from the pure flux-limited sample at different redshifts; for 
this reason, the quantitative measurements presented in 
this paper will all be based on quasi-volume-limited sam- 
ples, limited to an absolute magnitude Mb — —20+5 log h. 
Finally, it should also be mentioned that in adopting a uni- 
versal luminosity function we do not take into account the 
possible dependence of the luminosity function on mor- 
phological type and environment; this is, however, a sec- 
ond order effect in this work. 

The shot-noise error affecting the reconstructed field 
at different r is estimated by computing the square root 
of the variance 

T / „/IAr.l\ \ 2-, 1/2 

<r) = = 
Pg 



E 



S{n, M c )<f>(mi), C(z, to), (a, S) 



The amplitude of the shot noise increases as a function 
of redshift in a purely flux-limited survey. We deconvolve 
the signature of this nois e from the density maps by apply- 
ing the Wiener filter (cf. IPress et al , 1992 ; Strauss et al 



1992b) which provides the minimum variance reconstruc- 
tion of the smoothed density field, given the map of the 
noise and the a priori knowledge of the underlying power 
spectrum (e.g., Lahav et all Il994). For this we assume 



where brackets denote statistical averages and where 
P e (k) = (27r)~ 3 (|e 2 (k)|) is the power spectrum of the 
noise. Assuming ergodic conditions, this last quantity can 
be computed as P e (k) ~ (27r)~ 3 |e(k)| 2 . The calculation 
of ((5f.(k)) taking into account the form of the window 
function F and the peculiar VVDS survey geometry is 
presented in paper I. 



3.2. A cosmographical tour up to z = 1.5 

We have first applied our reconstruction technique to the 
global flux-limited VVDS sample to build a visual three- 
dimensional map of galaxy density fluctuations to z = 1.5 
which exploits the full information content of the survey. 
The / < 24 sample is characterised by an effective mean 
inter-particle separation of ((r) ~ 5.1 /i _1 Mpc ) in the red- 
shift range [0,1.5]. For comparison, this sampling is better 
(denser ) than the early CfAl survey ((r) ~ 5.5ft. _1 Mpc) 
used by lDavis k. Huchral ( 198lh to reconstruct the 3D den- 
sity field of the local Universe (i.e., out to ~ 80 /i _1 Mpc ). 
Also, at the median depth of our survey, i.e., in the redshift 
interval 0.7 < z < 0.8, the mean inter-particle separation 
is 4.4 h~ 1 Mpc, a value nearly equal to the 2dFGRS at its 
median depth. 



that the observed galaxy density field S g (r) , and the true Universe was 
(i.e. including all galaxies) underlying field ^r(r), both 
smoothed on the same scale, are related via 



The recovered galaxy density field is presented in Fig. 
1. Fluctuations have been smoothed on a scale R = 
2/i _1 Mpc. Only density contrasts with signal-to-noise ra- 
tio S/N > 2 are shown. 

A remarkable feature of this "geographical" exploration 
of the Universe at early cosmic epochs is the abundance of 
large-scale structures similar in density contrast and size 
(at least in one direction) to those observed by local sur- 
veys. In particular, it is tempting to identify qualitatively 
a few filament-like density enhancements bridging more 
condensed structures along the line of sight, although the 
survey transverse size is still too small to fully sample their 
extent. Nevertheless, it is interesting to notice that these 
apparently one-dimensional structures remain coherent 
over scales ~ 100/i _1 Mpc, separating low-density regions 
of similar size. Figs, one and two visually confirm that the 
familiar web pattern observed in the local Universe is not 
a present-day transient phase of the galaxy spatial organ- 
isation but it is already well-defined at ~ 1.5 when the 
30% its p resent age (e.g . 



1995 iGerke et all [200 



p resent age (e.g . , iLe Fevre et al 
I IScoville et all |2007j). This 7" 



6 g (r) = S T (r) + e(r), 



(5) 



where e(r) is the local contribution from shot noise (see 
Eq. ([4])). The Wiener filtered density field, in Fourier 
space, is 



5jr(k) = r(k)6 g (k) , 

where 

(4(k)> 



(<5 2 (k)) + (2 7 r)3p £ (k) 



(0) 



(7) 



plies that large-scale features of the galaxy distribution 
essentially reflects the long- wavelength modes of the initial 
power spectrum, in agreement with theoretical predictions 
of the CDM hierarchical scenario. Numerical simulations 
of large scale structure formation in fact show that the 
present-day web of filaments and walls is actually present 
when the universe was in embryonic form in the overden- 
sity pattern of the initial fluctuations, with subsequent 
linear and non-linear g ravitational dynamics just sharp- 
ening its features (e.g. 
l2005l i. 



Bond et al. . 1996; ISpringel et al 
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Fig. 1. The reconstructed density field for 0.4 < z < 1.4, as traced by the galaxy distribution in the VVDS-Deep 
redshift survey to / < 24. This figure preserves the correct aspect ratio between transverse and radial dimensions. The 
mean inter-galaxy separation of this sample at the typical depth of the VVDS (z = 0.75) is 4.6/i _1 Mpc, comparable 
to local redshift surveys as the 2dFGRS. The galaxy density distribution has been smoothed using a 3D Gaussian 
window of radius R = 2/i _1 Mpc and noise has been filtered away using a Wiener filtering technique (see Strauss & 
Willick 1995, Marinoni et al. 2005). Only fluctuations above a signal-to-noise threshold of 2 are sh own. The ac c uracy 



and robustness of th e reconstruction methods have been tested using realistic mock catalogues (jPollo et all 12005 



Marinoni et all 120051) . 
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The limited angular size of the survey is exemplified 
by a dense "wall" at z = 0.97 that stretches across the 
whole survey solid angle (0.7 x 0.7 deg) (see Fig. 2). 
This two-dimensional structure is coherent over more than 
~ 30/i _1 Mpc (comoving) in the transverse direction, is 
only ~ 10/i _1 Mpc thick along the line of sight, and has 
a mean overdensity 8 g = 2.4 ± 0.3. This makes it sim- 
ilar to the largest and rarest structures observed in the 
local Universe, such as the Shapley conce ntration (e.g. 
Scaramella et allll989UBardelli et alll2000h. By applying 



a Vor onoy-Dclaunay cluster finding code (M arinoni et al 



2002), we find 10 distinct groups in this structure, with 
between 5 and 12 galaxy members each (down to the lim- 
iting magnitude 1=24), for a total of 164 galaxies. If one 
considers the evolution of mass fluctuations in the stan- 
dard ACDM model, the probability of finding a struc- 
ture with similar mass overdensity at such early times 
(0.9 < z < 1) would be nearly 4 times smaller than to- 
day: one such mass fluctuation would be expected in a 
volume of <~ 3 ■ 10 6 /i~ 3 Mpc 3 , i.e., nearly 5 times larger 
than our surveyed volume up to z ~ 1. In fact, as we shall 
describe in section 3.3, finding such a galaxy overdensity 
is not so unusual: it is clear evidence that the biasing be- 
tween galaxies and matter at these epochs is higher than 
today. This makes fluctuations in the galaxy distribution 
to be highly enhanced with respect to those in the mass. 



3.3. Evolution of the PDFs of galaxy fluctuations in 
the VVDS 



Several approaches may be used to characterise in a quan- 
titative way the distribution of galaxy fluctuations S g 
shown in Fig 1. A complete specification of the overdensity 
field may be gi ven by the full set of gal axy N-point correla- 
tion functions ( Davis fe Peebles! . 19771 ). This approach has 
been explored and routinely applied over the past decade 
as better and deeper redshift surveys have become avail- 
able. An alternative description may instead be given in 
terms of the probability distribution function of a ran- 
dom field. By definition, the PDF of cosmological density 
fluctuations describes the probability of having a fluctu- 
ation in the range (5,6 + dS), within a spherical region 
of characteristic radius R randomly located in the survey 
volume. In principle, it encodes all the information con- 
tained within the full hierarchy of correlation functions, 
and provides insights about the time evolution of density 
fluctuations. This definition can be applied either to the 
distribution of galaxies, characterizing their number den- 
sity fluctuations, or to the dark-matter dominated mass 
distribution. For the latter case, the expected shape of 
the PDF can be predicted as a function of redshift given 
a cosmological model, at least for large-scale fluctuations; 
this can be done either analytically (see below) or using 
numerical simulations. 

On the observational side, in the case of surveys of the 
local Universe this fundamental statistics has been often 
overlooked (but see Marinoni & Hudson 2002, Ostriker 



et al. 2003). On the other hand, only recently deep red- 
shift surveys have reached sufficient volumes to allow these 
measurements to be extended back in time. In paper I, we 
have discussed and tested in detail the methodology to 
estimate the PDF from this kind of samples. In particu- 
lar we have checked the robustness of the reconstruction 
against the specific VVDS-02h survey selection function, 
shot-noise errors and other observational biases. We used 
fully realistic mock samples of the VVDS-02h survey data 
and showed that, once the smoothing scale R is larger 
than the mean inter-galaxy separation, the overall shape 
of the reconstructed PDF is an unbiased realisation of 
the complete parent galaxy population. In particular, we 
showed that for redshifts up to z = 1.5 the VVDS-02h sky 
coverage and sampling rate are sufficient for obtaining a 
reliable reconstruction of the PDF shape (in both low- and 
high-density regions) on scales R > 8/i _1 Mpc. Clearly, the 
degree with which the PDF measured from this sample is 
a fair representation of the "universal" PDF up to z=1.5 
is a separate, yet critical question. A difference is natu- 
rally expected due to fluctuations on scales larger than 
the volume probed ("cosmic variance"). An estimate of 
this effect is actually included in our error bars, as these 
were drawn from the scatter among our set of VVDS mock 
samples. 

We have therefore applied the estimator of Eq. [2] and 
the full dc- noising technique described in § 13.11 to the two 
luminosity-limited sub-samples of our survey described in 
§ [2j reconstructing the PDF of galaxy fluctuations in top- 
hat spheres of radius R — 10/i~ 1 Mpc at two different 
epochs (0.7 < z < 1.1 and 1.1 < z < 1.5). The typical 
luminosity of the galaxies selected in these two intervals 
(Mb < — 20 + 5 log h) corresponds to a median luminosity 
Lb — 2L|j at z ~ (i.e. th e same median lum inosity of 
the whole 2dFGRS sample (| Verde et all l2002h V As dis- 
cussed previously, the use of luminosity-selected samples 
virtually eliminates distance-dependent shot-noise contri- 
butions (clearly neglecting the residual evolution within 
the two redshift bins, which is well within the errors). 
The results are shown in Fig. 3. The measured PDF's in 
Fig. 3 show several interesting features. First, as times 
passes (redshift decreases) the maximum of the PDF shifts 
to smaller <5-values; secondly, the low-density tail is en- 
hanced, with more low-density regions appearing at lower 
redshifts. Quantitatively, this implies in particular that 
the probability of having an under-dense (S g < 0) region 
of radius R = 10/i _1 Mpc at 0.7 < z < 1.1 is nearly 10% 
larger than at earlier times (1.1 < z < 1.5). 

3.4. The expected PDF of mass fluctuations in the 
mildly non-linear regime 

The shape of the PDF of the galaxy overdensities is 
strongly dependent on the non-linear effects implicit 
both in the gravitational growth and in the physi- 
cal mechanisms respo nsible for galaxy formation (e.g. 
Watts fe TavlorL 120001 ) . Initial density fluctuations are 
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Fig. 2. Density distribution and properties of a large-scale planar structure at z = 0.97, that completely fills the 
VVDS-02h field-of-view. 



normally assumed to have a Gaussian PDF; this is then 
modified by the action of gravity and, in the case of the 
galaxy field, by the way galaxies trace the underlying mass 
(biasing scheme). If galaxies were faithful and unbiased 
tracers of the underlying mass, the peak shift and the de- 
velopment of a low-density tail we observe in Fig. 3 could 
be naturally interpreted as the key signature of dynamical 
evolution purely driven by gravity. In fact, gravitational 
growth in an expanding Universe makes low density re- 
gions propagate outwards and become more common as 
time goes by, while at the same time the high-density tail 
increases. 



If this interpretation is correct, we expect the PDF of 
galaxy overdensities to coincide with the PDF of mass fluc- 
tuations in each redshift range, once they are normalised 
to the observed clustering at z ~ 0, where we know that 
L ~ 2L* galaxies trace the mass (jVerde et all 120021) . Let 



us verify whether this is the case by first summarising the 
main formalism to compute the PDF of mass fluctuations 
in a given cosmological scenario. 

In hierarchical models, it is well established from nu- 
merical simulations that when structure growth reaches 
the nonlinear regime on a scale R, the PDF of mass den- 
sity contrasts in comoving space is well described by a log- 



normal distribution rtColes fc Jonesl 1991 L Kofman et al. , 



1994HTaylor fc WattkbOOOtlKayo. Taruva k Sutol . l200lh . 



f R (S) = 



(2tt^|)- 1 / 2 - _ f [ln(l + S)+ c4/2] 2 



1 + 8 



exp 



{ 



2c4 



} • (8) 



This is fully characterised by a single parameter u)r, re- 
lated to the variance of the (5-ficld on a scale R as 



u% = Hi + (s 2 ) R ] 



(9) 



At high redshifts, the variance a 2 R = (S 2 )r over suf- 
ficiently large scales R (those explored in this paper) is 
given in the linear theory approximation by 



<tr{z) = <?r(0)D{z) 



(10) 



where D(z) is the linear growth factor of density fluctua- 
tions (normalised to unity at z = 0), 



D(z)=exp-[/ f(z)dln(l + z) 
'-Jo 



(11) 



In the standard ACDM cosmological model, the expres- 
sion for the logarithmic derivative of the growth factor, 
f = dlogD/dloga (with a = (1 + z)^ 1 ), can be approxi- 
mated to excellent accuracy as 

f(z) ~ Q&(x) 
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where 7 ~ 0.55 (Wang & Steinhardt 1998, Linder 2005) 
and 

n (z )-np ( 1 + ^) 3 
E 2 (z) = n° m (i + zf + n A . 

The lognormal approximation formally describes the 
distribution of matter fluctuations computed in real co- 
moving coordinates. On the contrary, the PDF of galax- 
ies is observationally derived in redshift space, where 
its s hape is distorted by the effects of peculia r motions 



(e.g. iMarinoni et all IT9981 IGuzzo et all l2008h . In order 



to map properly the mass overdensities into galaxy over- 
densities the mass and galaxy PDFs must be computed 
in a common reference f rame. It has been shown by 
Sigad. Branchini fc Dekell (2000) that an optimal strat- 
egy to derive galaxy biasing is to compare both mass and 
galaxy density fields directly in redshift space. Implicit in 
this approach is the assumption that mass and galaxies 
are statistically affected in the same way by gravitational 
perturbations, i.e. that there is no velocity bias in the mo- 
tion of the two components. 

The relation between the variances measured in real 
and redshift comoving space is 



OflO) =P{z)a R {z) 



(12) 



where p(z) is a redshift-dependent correcting factor which 
takes into account the average contribution of th e linear 
redsh ift distortions induced by peculiar velocities (jKaiser , 
Il987|). Its expression, in the high redshift regime, is given 



by ([Hamilton! . Il998t IMarinoni et all 120051 ) 



^)=[l + ^(z) + if 2 (z)] 1/2 



(13) 



We have used this formalism to compute the PDF ex- 
pected for mass fluctuations in the redshift ranges ex- 
plored using our galaxy samples. This is given, for the 
two ranges, by the curves in the top panel of Fig. 3. The 
evident discrepancy between the galaxy and mass PDF's 
indicates that the observed evolution cannot be only the 
product of gravitational growth (in the adopted cosmo- 
logical model), but that a time-evolving bias between the 
galaxy and mass density fields is needed: at high redshifts 
and on large scales galaxy overdensities trace in a more 
biased way the underlying pattern of dark matter fluctu- 
ations. In the following section we shall summarise our 
current knowledge on the properties and evolution of the 
biasing function and show how the presence of a non-linear 
bias is a necessary ingredient to theoretically understand 
the evolution of the PDF in Fig. 3. This will in particular 
provide us with the necessary background to interpret the 
evolution of the low-order moments of the PDF at differ- 
ent redshifts, which is our aim in this paper. 




■1.5 -1.0 -0.5 0.0 0.5 1.0 
log(y) 




0.00 

-1.5 -1.0 -0.5 0.0 
log(y) 

Fig. 3. The PDF of galaxy fluctuations (in units y = 1+S) 
for VVDS galaxies with M B < -20 + 5 log h from the 
VVDS within two independent volumes, corresponding to 
different cosmic epochs: 0.7 < z < 1.1 (blue shaded his- 
togram), and 1.1 < z < 1.5 (green shaded histogram). 
The galaxy PDF has been reconstructed using a Top- 
Hat smoothing window of comoving size R = 10/i _1 Mpc. 
The histograms actually correspond to the distribution 
function G(y) = ln(10)y<?(y) because the binning is done 
in log(y). The two observed histograms have been repro- 
duced in both the upper and lower panels. They are com- 
pared to the theoretical predictions for the PDF of, respec- 
tively, mass fluctuations (top, from Eq. [5]) and of galaxy 
fluctuations as inferred from Eq. [TJ] using the non-linear 
biasing function measured from the VVDS (bottom). The 
blue and red lines correspond to the higher- and lower- 
redshift samples respectively. 



3.5. Evolution and non-linearity of biasing 

Biasing lies at the heart of all interpretations of large-scale 
structure models. Structure formation theories predict the 
distribution of mass; thus, the role of biasing is pivotal in 
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mapping the observed light distribution back into the the- 
oretical model. In our case we need to disentangle the im- 
print of biasing from that of pure gravity in the evolution 
of the galaxy PDF. 

In paper I we inferred the biasing relation S g — S g (5) = 
b(S)8 between mass and galaxy overdensities from their re- 
spective probability distribution functions f(8) and g(S g ). 
Assuming a one-to-one mapping between mass and galaxy 
over density fields, conservation of probabilities impl ies 



over aensity nelds, con servation ot probabil ities impl i 
(e.g. ISigad. Branchini fc DekeJ . l200(j IWild et all l2005h 



g(S g )dS g = f(5)dd 



(14) 



This approach implies the assumption of a cosmologi- 
cal model (the standard ACDM model in our case) from 
which to compute f(5), the mass PDF. The advantage 
over other methods is that we can explore the functional 
form of the relationship S g = b{z, 6, R)S over a wide range 
in mass density contrasts, redshift intervals and smoothing 
scales R without imposing any a-priori parametric func- 
tional form for the biasing function. Note that, by defi- 
nition, this scheme is ineffective in capturing information 
about possible stochastic properties of the biasing func- 
tion. 

The numerical solution 5 g = 5 g (5) of Eq. [14] maps the 
mass PDF (solid lines in the top panel of Fig. 3) into 
the galaxy PDF (solid lines in the bottom panel of Fig. 
3) and can be analytically approxi mated using a Taylor 
expansion (|Frv fc Gaztahaga , 1993) 



(15) 



fe=0 



where the coefficients bi depend on redshift. We consider 
this power series only to second order, and fit the nu- 
merical solution for the biasing function leaving bo as a 
free parameter. Avoiding setting bo as an integral con- 
straint (< S g >= 0) allows us to account for possible (un- 
modelled) contributions from higher order moments of the 
expansion. This approach has the advantage of minimis- 
ing biases in the estimates of the lower moments of the 
expansion, specifically b\ and &2- 

The key result from Paper I has been to show that 
galaxy biasing is poorly described in terms of a sin- 
gle scalar and better characterised by a more sophis- 
ticated representation. Specifically, always considering a 
scale R — 10/i~ 1 Mpc, the ratio between the quadratic 
and linear bias terms has been evaluated in four different 
high redshift intervals (see Table 2 of Paper I). When av- 
eraged over the full redshift baseline 0.7 < z < 1.5, this 
ratio turns out to be 



(16) 



■f- ) = -0.19 ±0.04 



i.e. different from zero at more than 4cr confidence level. 
This means that - at least over the redshift range and 
scales considered here - the level of biasing depends on 
the underlying value of the mass density field. In other 
words, the way galaxies are distributed in space depends 



in a non-linear manner on the local amplitude of dark 
matter fluctuations. 

The measurement of a non-linear term in the biasing 
relation is fully consistent with a parallel analysis of the 
hierarchical scaling of the N-point correlation functions 
in the same VVDS sample (Cappi et al. 2008). These re- 
sults confirm a gener ic prediction of hierarc hical models 



of galaxy formation ( Somerville et al. . l200ll e.g., ). It is 



relevant to compare them to estimates of the bias func- 
tion at the current epoch. Early works indirectly suggested 
that also at z ~ the biasing function should have a 
non-negligible non- linear component. Comparing the two- 
point correlation functions and the normalized skewness 
of SSRS2 galaxies, Benoist et al. (1999) showed that the 
relative bias between galaxies with different luminosity is 
non-linear, which indirectly indicates that (at least for lu- 
minous galaxies) the bias with respect to the dark matter 
must be non-linear as well. A similar analysis was per- 
formed by Baugh et al. (2004) and Croton et al. (2004) on 
the 2dFGRS, finding results consistent with Benoist et al.; 
finally, Feldman et al. (2001) and Gaztanaga et al. (2005) 
directly measured the three-point correlation function (in 
both Fourier and real space) for the IRAS and 2dFGRS 
samples respectively and found evidence for 62 < 0. 

On the other hand, these results seem to be incon- 
sistent with another a nalysis of the 2dFG RS performed 
using the bi-spectrum ( Verde et al. . 1200*21 ). Hikage et al. 
2005 by analysing the SDSS galaxies with a Fourier- 
phase technique also conclude that the bias in this sur- 
vey is essentially linear. If one ignored the other inde- 
pendent analysis of the 2dFGRS, it could be speculated 
that the large-scale non-linear term that we detect at 
z > 0.7 is suppressed as a function of cosmic time; this 
is however not supported by the r esults of numerical ex- 



periments (jSomerville et all 120011 ). Interestingly, we note 



that when compa red to the local non-linear measurements 
(&2/ &1 ~ —0.35) (jFeldman et all l200ll ; iGaztanaga et all 
2005), our estimate seems to suggest that the amplitude 



of the quadratic term 62/&1 decreases (in absolute terms) 
as a function of redshift, a results in qualitative agreement 
with indications from simulations. It seems therefore more 
likely that the difference in the reconstruction methods 
used (with different sensitivity to higher order terms in 
Eq. [TBf is a better explanation of the discrepant results 
at z ~ 0. We will show in the next section (§4) how the 
self-consistency of the evolution of the variance and skew- 
ness of galaxy counts with redshift, indeed requires the 
presence of a non- linear biasing component. 

The information contained in the non-linear function 
of Eq. Qjj] can be compressed into a single scalar term 
that can be used to interpret the evolution of two-point 
statistics (correlation function) as well as the variance of 
the galaxy density field (see §4.2). Since, by definition, 
(b(5)5) = 0, the most interesting linear bias estimators 
are associated to the second order moments of the PDFs, 
i.e., the variance (S g ) and the covariance (S g 5). Following 
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the prescriptions of Dekel &: Lahav ( 1999h . we character- 
ize the biasing function as follows 



(b 2 {5) 5 2 



(17) 



where h\ is an estimator of the linear biasing parameter 
defined, in terms of the two-point correlation function, as 
i g = b\i. We evaluate Eq. [T7] using Eq. [15] with parame- 
ters bi(z) estimated locally by Verde et al. (2002) and in 
the rcdshift range 0.4 < z < 1.5 by Marinoni et al. (2005). 
The best fitting phcnomcnological model describing the 
redshift scaling of the linear biasing parameter for a vol- 
ume limited population of "normal" galaxies with median 
luminosity L~2L*(z = 0)is 



b L (z) = l + (0.03±0.01)(l + z) 



3.3±0.6 



While today ~ 
mass distribution 



2L* galaxies tr ace the und e rlyin 



on large scales (jLahav et all 1200 



Verde et all 120021 ; iGaztanaga et al. . 2005 ). in the past the 
two fields were progressively dissimilar and the relative 
biasing systematically higher. In Paper I we showed how 
this observed redshift trend compares to different theoret- 
ical models for biasing evolution, i.e. a "galaxy conserv- 
ing" model (Fry et al. 1996), a "halo merging" model (Mo 
& White 1996) and a "star forming" model (Tegmark & 
Peebles 1998). 

4. Testing gravitational instability with the 
low-order moments of the PDF 

Having decoupled biasing effects from the purely gravita- 
tional evolution of the galaxy PDF we have now all the 
ingredients to use this latter quantity to test the consis- 
tency of some general predictions of the GIP. The evo- 
lution of the low-order statistical moments of the galaxy 
PDF, specifically its second and third moments can be 
compared, on large scales with analytical predictions of 
linear and second order perturbation theory respectively. 

4.1. Estimating the moments from redshift survey data 

Following standard conventions, we define the second- and 
third-order moments, on a scale R, of a continuous, zero- 
mean overdensity field as 



and 



<*2>J 



Sg9R(Sg)dS g . 



Sg9R(Sg)dS g . 



(18) 



(19) 



Note that the moments cannot be estimated as ensem- 
ble averages over the reconstructed PDF. In fact, this last 
quantity has been reconstructed using the Wiener filter- 
ing technique. This minimises the shot noise contribution 
(§ 13. ip but gives a biased estimate of the density field mo- 
ments (via Eq. [18] and fl9|) as it requires an input power 



spectrum (and therefore assuming a second moment). A 
standard practical way to estimate moments is to ran- 
domly throw spherical cells down within the galaxy distri- 
bution and reconstruct the count probability distribution 
function Pk = rik/N (where rik is the number of cells that 
contain k galaxies out of a total number of cells N. The 
moments are then estimated as 



{S p g ) = N- p J2Pk{(k~ N) P ) 



(20) 



k=0 



where N = J2T=o kp k- 

The quantities we are interested in are the cumulants 
(S p ) c of the one-point density PDF. For a density field 
smoothed with a top-hat window, the p-order cumulant 



1 



£ p (r ll r 2 ...r p )d 3 r 1 d 3 r 2 ...d 3 r p 



(21) 



is the average of the N-point reduced correlation func- 
tion over the corresponding cell of volume vr (from now 
on we will only consider the scale R = 10/i _1 Mpc and we 
will drop the suffix R, unless we need to emphasize it). 
This is defined as the connected part of the N-point cor- 
relation function (Sg(ri)5 g (r2)...Sg(r p )) in such a way that 
for p > 2 £ p = for a Gaussian field. Since the galaxy dis- 
tribution is a discrete process (Eq. [2] is a sum over Dirac 
delta functions) and since, by definition, the density con- 
trast has a zero mean, the connection between low-order 
cumulants and moments is given by 



1 



(22) 



(23) 



These relations accounts for discr eteness effects using 
the Poisson shot-noise model (e.g., IPeeblesL Il98fi IfTvL 
1985). Possible biases introduced by this technique are 
discussed by Hui & Gaztanaga 1999, while an alternative 
approach is detailed by Kim & Strauss (1998). 

Finally, it is necessary to devise a strategy to compen- 
sate for the fact that a cell will sample regions that have 
varying angular and spectroscopic completeness and which 
may even span the survey boundary. For this reason the 
galaxy counts are scaled up in proportion to the degree 
of incompleteness in the cell. This is done by weighting 
galaxy counts using the selection functions $(m), £(r, to), 
and 'J' (a, S) defined in ij l3.ll Additionally, although our re- 
construction scheme accounts for the non-uniform VVDS 
angular coverage, we further restrict ourselves to counts 
in spheres having at least 70% of their volume in the 
denser 4-pass region, in order to avoid possible edge ef- 
fects. We remark that moments are estimated from virtu- 
ally volume-limited samples, as defined in §[2] As a conse- 
quence, the radial selection function is constant and any 
variations in the density of galaxies are due only to large- 
scale structure. 
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4.2. The evolution of the rms and skewness of galaxy 
fluctuations 

Since in perturbation theory higher order cumulants are 
predicted to be a function of the variance, it is useful, in 
the following, to define the normalized skewness 



s 3 = (5l) c /a 2 



(24) 



where the shot-noise corrected variance a 2 is given by 
Eq. [22] Fig. 4 shows the evolution of the rms fluctuation 
and the normalized skewness on a scale R = lO/i^Mpc, 
as measured from the VVDS volume-limited sub-samples. 
Errors have been computed using the 50 fully-realistic 
mock catalogs of VVDS-Deep discussed in Polio et al. 
(2005). This allows us to include an estimate of the con- 
tribution of cosmic variance, which represents the most 
significant term in our error budget. 

The top panel of Fig. 4 shows that the square-root 
of the variance, which measures the r.m.s. amplitude of 
fluctuations in galaxy counts, is with good approxima- 
tion constant over the full redshift baseline investigated: 
in redshift space, the mean value of a g for our volume- 
limited galaxy samples is 0.78 ± 0.09 for 0.7 < z < 1.5. 
A similar, nearly constant value is also consistent with 
the value estimated at z ~ 0.15 f rom the 2dF galaxy red- 



shift survey (jCroton et al.l . 120041) that is also reported in 



same figure. This means that over nearly 2/3 of the age 
of the Universe the observed fluctuations in the galaxy 
distribution look almost as frozen, despite the underlying 
gravitational growth of mass fluctuations. This quantifies 
the visual impression we had from Fig. 1, that the distri- 
bution of galaxies is as inhomogeneous at z ~ 1 as it is 
today. 

The third moment, which measures asymmetries be- 
tween under- and over-dense regions, indicates that the 
galaxy density field was non-Gaussian on large scales (10 
h~ 1 Mpc ) even at these remote epochs (~ 4er detection). 
In particular we find indication for an increase of the nor- 
malised skewness with cosmic time, when comparing the 
VVDS values to the local measurement by 2dFGRS. 

Using the measured bias evolution, we can translate 
the specific predictions of the GIP for the variance and 
skewness of the matter density field into the corresponding 
observed quantities. Using linear perturbation theory, the 
scaling of the rms of number density fluctuations is 



<j g (z) ~ b L (z)D(z)p(z)a(0) 



(25) 



In a Universe in which primordial density fluctuations 
were Gaussian, the non-linear nature of gravitational dy- 
namics leads to the emergence of a non-trivial skewness of 
the local density PDF. 

Within the framework of gravitational perturbation 
theory, the first non-vanishing term describing the evolu- 
tion of the skewness of a top-hat filtered, initially Gaussian 
matter density field corresponds to second-order. 

According to non-linear, second-order perturbation 
theory predictions, the skewness of the mass distribu- 



R=10h"'Mpc 




Fig. 4. Evolution of the r.m.s deviation (top) and skew- 
ness (bottom) of the PDF of galaxy fluctuations on a 
scale R = 10/i _1 Mpc. The filled squares correspond to 
two volume-limited samples from the VVDS with Mb < 
—20 -I- 51og/i covering the redshift intervals indicated by 
the shaded regions. Triangles correspond to the 2d FGRS 



measurements at z ~ 0.15 ([Croton et al.l . 120041 ). from 



a sample including similarly bright galaxies. Error bars 
give 68% confidence errors, and, in the case of VVDS 
measurements, include the contribution from cosmic vari- 
ance. The dashed lines in both panels show the theo- 
retical predictions for the evolution of the variance (Eq. 
251) and skewness (Eq. [27|) inferred using VVDS mea- 
surement of biasing. Predictions for the skewness (based 
on the (b\{z),b2{z)) measurements in the redshift range 
0.7 < z < 1.5 quoted in Table 2 of Paper I) have been 
extrapolated to z ~ using the local (2dFGRS) biasing 
measurements of Verde et al. (2002) (linear bias, dotted 
line) and of Gaztahaga et al. (2005) (quadratic bias with 
= —0.34, dot-dashed line). 

tion is approximately independent of time, scale, den- 
sity, or geometry of the cosmological model. Assuming 
that its evolution only depends on the hypothe- 
sis that the initial fluctuations are small and quasi- 
Gaussian and that they grow via gravitational clus- 
tering o ne der ives that, in redshift-distorted space 



(Peebles, 1980; Juszkicwicz, Bouchct . fc Colombil . 11993 



BernardeauL Il994t Ifflvon et all Il995^ 



S, 



35.2 



1.15(ra + 3) 



(26) 



where n is the effective slope of the power spectrum on the 
scales of interest [i.e. in our case, sinc e R = 10fe~ 1 Mpc , 



n is approximately given by -1.2 (e.g. iBernardeau et al.l . 
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20021) ]. Subsituting the evolution of bias in second order 
approxim ation, the evolution of th e observed skewness is 
given by (|Frv &; Gaztaiiagal Il993l ) 



5 3,g 



S, 



(27) 



The curves in both panels of Fig. 4 show that equations 
(|25|) and (|27jl reproduce extremely well the evolution of 
variance and skewness observed within the VVDS. 

The mass PDF is a one-parameter family of curves 
completely specified once the linear evolution model for 
the mass variance (S 2 ) is supplied. This implies that our 
non-linear biasing estimate is fully independent from pre- 
dictions of higher-order perturbation theory. On the con- 
trary, non-linear biasing at z — is inferred by directly 
matching 3-point galaxy statistics with the correspond- 
ing mass statistics derived from weakly non-linear per- 
turbation theory (e.g. Verde et al. 2000, Gaztahaga et 
al. 2005). As a consequence, the agreement we find be- 
tween predicted and observed third-order moments is not 
a straightforward consequence of the method used to de- 
rive the biasing function. These results provide an indi- 
cation of the consistency, at z = 1, of some constitutive 
elements of the standard picture of gravitational instabil- 
ity from Gaussian initial conditions. 

Concerning the local measurements from 2dFGRS, the 
predicted scaling for the skewness continues to show very 
good agreement if the local, non-linear measurement of 
Gaztahaga et al. (2005) is considered. The value of £3,3, 
however, cannot be consistent with GIP predictions if in 
the local universe the simple linear biasing measurement 
of Verde et al. (2002) (i.e. b 2 = 0) is adopted. 

4.3. Effect of galaxy luminosity evolution 

In the above comparison of galaxy samples at three differ- 
ent epochs, we have so far neglected an important point. 
Galaxy luminosity evolves significantly between z ~ and 
z ~ 1, with a mean brightening of at least 1 magnitude 
for an average spectral type (Ilbert et al. 2005). Thus, 
the contribution to the clustering signal at progressively 
earlier epochs may not be be due to the progenitors of 
the galaxies that are sampled at later times in the same 
luminosity interval. Luminosity evolution between 2 = 1 
and z = 1.5 is more uncertain but certainly smaller due 
to the shorter time interval. The brightening between the 
two VVDS sub-samples considered is expected not to be 
very significant for these very luminous objects. 

To compare galaxies at z = 1 and z ~ one should in 
principle confront our high-redshift results with those of 
local galaxies about one magnitude fainter. According to 
2dFGRS results, this means shifting the local estimates of 
variance to a slightly smaller value (Norberg et al. 2002) 
and leaving the skewness measurement unchanged within 
the quoted error bars (Croton et al. 2004). These changes 
would only reinforce our conclusions about the evolution 
of the low order moments of the PDF of galaxy fluctua- 
tions. In particular, since a fainter sample has a smaller 



bias threshold, the locally measured skewness would make 
the discrepancy with GIP predictions for a simple linear 
biasing model even stronger. 



5. Conclusions 

The results presented in this paper provide the hrst direct 
evidence at z ~ 1 for the consistency of the GIP hypothe- 
sis as described in the framework of general relativity. The 
standard theory of structure formation via gravitational 
instability successfully explains the present day statistics 
(e.g. Tegmark et al. 2006) and dynamics (e.g. Peacock et 
al. 2001) of large scale structures. We have shown that ob- 
servations arc fully consistent with these predictions over 
the entire redshift baseline < z < 1.5 once the biasing 
between the galaxy and matter distribution is properly 
described. In Paper I we showed that it is necessary to 
include a small (10%) yet crucial non-linear component to 
accurately account for the observed probability distribu- 
tion function of galaxy overdensities. Here we have shown 
that this component is also necessary as to understand 
the observed evolution of the low-order moments of the 
galaxy overdensity field. 

More specifically, our analysis of the 3D density fluc- 
tuation field traced by a volume-limited sample of VVDS 
galaxies (with Mb < —20 + 5 log ft.) at different epochs 
unambiguously reveals the time-dependent effects of grav- 
itational evolution: 

a) underdense regions progressively occupy a larger 
volume fraction as a function of cosmic time, as expected 
from gravitational growth in an expanding background; 

b) the second moment of the field traced by this 
"normal" population of galaxies (with median luminos- 
ity ~ 2L*) is statistically consistent with the local (z ~ 0) 
estimate for similarly luminous galaxies, i.e. it is approx- 
imately constant over the full redshift baseline < z < 
1.5. This implies that the apparent inhomogeneity in the 
galaxy distribution remains similar, i.e. galactic fluctu- 
ations have alm ost frozen over nearly 2/3 of the age 



of the universe (iGiavalisco et al. . 19981 : Coil et al. . 2004 : 
23). 1 



iPollo et all 12005). We have shown that this is readily ex- 
plained by the combination of the gravitational growth of 
mass fluctuations with the evolution of the bias between 
galaxies and mass. These two factors almost cancel each 
other out; 

c) there are some hints that the skewness increases 
with cosmic time, its value at z ~ 1.5 being nearly 
2cr times lower than that measured locally by the 2dF- 
GRS for similarly luminous galaxies. In particular, the 
measured value of the skewness at z ~ 1.5 (on scales 
R = 10/i _1 Mpc ) indicates that galaxy fluctuations are 
strongly non-Gaussian (~ 4cr detection) even at such an 
early epoch (see Cappi et al. 2008 for a different approach 
which arrives at similar conclusions); 

d) remarkably, once VVDS measurements of non-linear 
biasing are included, both these trends are consistent with 
predictions of linear and second-order perturbation the- 
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ory for the evolution of gravitational perturbations as de- 
scribed within the framework of general relativity; 

e) we have shown that the values of the skewness 
we measure at high redshift are difficult to reconcile 
with the 2dFGRS measurements if local biasing is linear 
( Verde et al. . 20021) . A fully coherent gravitational picture 
emerges, however, over the whole baseline < z < 1.5 if 
the non-linearity of the local biasing function is taken into 
account, at the level estimated by Feldman et al. 2001 and 
Gaztahaga et al. 2005. Compared to these local measure- 
ments our results seem to suggest that the amplitude of 
the quadratic term 62/^1 is a decreasing function of red- 
shift at least up to z ~ 1.5. 
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